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Abstract 

HEPGEN is a generator of Monte Carlo events, which is dedicated to studies of hard exclusive 
single photon or meson production processes at the COMPASS experiment kinematic domain. In 
addition HEPGEN allows to generate also single photon or meson production accompanied by 
the diffractive dissociation of the nucleon, which is one of the main sources of background in 
analyses of exclusive processes. 
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1 Introduction 



HEPGEN (Hard Exclusive Production GENerator) is a dedicated Monte Carlo generator of events, 
which is used for studies of hard exclusive leptoproduction processes at the COMPASS experiment 
kinematic domain. The studied reactions comprise both single photon production via Deeply Virtual 
Compton Scattering (DVCS) and Bethe-Heitler (BH) process as well as hard exclusive production of 
various mesons (HEMP). In addition to exclusive processes, HEPGEN allows also to generate events 
of single photon or meson production accompanied by the diffractive dissociation of the nucleon, 
which is one of the main sources of background in analyses of exclusive processes. 

Presently there are four exclusive processes implemented in the generator: 

(a) single photon production (DVCS + BH), 

(b) exclusive 7r° production, 

(c) exclusive p° production, 

(d) exclusive p + production. 

The implementation of two other exclusive processes, 4> and J/ip production, is under development. 

2 Modelling of exclusive processes 

The main features of modelling of the specific processes are shortly described in this section. 

(a) The DVCS amplitude is described following the model of Frankfurt, Freund and Strikman 
(FFS) JT||2l with modifications required for COMPASS [3J. The FFS model was proposed to describe 
DVCS data at very small x Bj (~ 1CT 3 ) measured at HERA and it relates the imaginary part of forward 
DVCS amplitude to the structure function F%. In fact the model accounts only for the contribution 
from the GPD H, which is dominant for the unpolarised DVCS cross section at small xbj- The 
modifications and extensions introduced for COMPASS include in particular: an additional option 
for a non-factorizable (Reggeized) t-dependence of DVCS amplitude, calculation of the real part of 
the amplitude using derivative dispersion relations, and a proper choice of F 2 parameterisation JH 
suitable for the COMPASS kinematic range. 

The Bethe-Heitler amplitude is described using formalism of Ref. [0 with the replacement of the 
approximate expressions for the lepton propagators P 1 and P 2 by the exact formulae [6]. For the 
Dirac and Pauli electromagnetic form factors the dipole formulae are used, which are adequate for the 
investigated region of small t ( < 1 (GeV/c) 2 ). 

Similarly, the DVCS-BH interference term is described following the formalism of Ref. [5] with 
the DVCS amplitude as explained above. 

(b) Exclusive n° production is simulated using cross sections for longitudinal (a^) and transverse 
(cry) virtual photons according to the GPD model of Goloskokov-Kroll 0. In this model the contribu- 
tion from the transversity GPD (with helicity flip of exchanged quark) is also included. A reasonable 
agreement is achieved with low energy data on exclusive n° production from JLAB. The generator 
calculates weights of generated events using a grid of a L and a T cross sections as functions of W, Q 2 
and t'. 

(c) Exclusive p° production is simulated using the shapes of kinematic distributions as a function 
of u, Q 2 and p^, which were measured by the NMC [8 J for the deuterium target in a kinematic range 
similar to that of COMPASS. The cross section is rescaled using the value from the GPD model of 
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Ref. [HI at a selected kinematic point W and Qq. The ratio R = a L /a T is parameterised using all 
available measurements from the range 5 < W < 300 GeV. 

(d) Similarly as for n° meson, exclusive p + production is simulated using predictions from the 
GPD model for this reaction by Goloskokov and Kroll [10]. The grid of a L and a T as a function of 
W, Q 2 and t' is used to calculate cross section for each generated event. 



3 Diffractive dissociation of the target nucleon 

HEPGEN offers an option to generate also events of leptoproduction of single photon or meson acc- 
ccompanied by diffractive dissociation of the target nucleon into several particles, mostly a nucleon 
and pions. 

As no experimental data exists for such processes at COMPASS kinematic range, a factorisation 
ansatz is used. We assume that the cross section for such a process factorises into a part related to 
the virtual-photon vertex and that for the nucleon vertex. The former is parameterised using available 
data or models for exclusive leptoproduction. For the latter one uses parametrisation of experimental 
data on the proton-proton and proton-deuteron scattering at fixed target Fermilab experiments, which 
were carried out with the proton beam energies in the range 100 - 300 GeV ifiTll . 

The following results from these experiments are used in HEPGEN. 

• The energy dependent ratio of the cross section for single diffractive dissociation, pp(d) — > 
Xp(d), to that for elastic scattering, pp(d) — > pp(d). 

• The shape of the cross section for single diffractive dissociation as a function of Mx, where Mx 
is the invariant mass of the diffractively produced system X. 

• Charged and neutral hadron multiplicities for the system X. 



4 HEPGEN output 

The HEPGEN generator runs in a standalone mode producing a file with the output data. This file 
can be read-in by COMGEANT, the GEANT-based Monte Carlo package to simulate responses of 
COMPASS detectors. The format of the HEPGEN output data is conformal with the structure of 
LEPTO/JETSET output. In particular, it contains kinematic information on all initial, intermediate 
and final particles generated in each event. Each event is assigned a weight, which is proportional to 
the cross section corresponding to the values of kinematic variables for that event. 

The information from the generator is propagated in its original format through the apparatus sim- 
ulation and event reconstruction stages, up to the final mini-DST files, which are used for the physics 
analyses. Thus this information could be accesses either immediately after HEPGEN run, or at the 
end of the simulation and generation chain. 

A detailed information on the installation of the HEPGEN package, required input files and steer- 
ing parameters data, as well as on how to access the generator output information is given in the 
following sections and can be found in Ref. |[T2l . 
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5 Technicalities 



In this section described are installation and usage of HEPGEN generator. 



5.1 Installation 

The latest version of HEPGEN source code can be download from the web page given as by Ref. lfl~2ll. 
The source code files are compressed with GNU zip compression utility. Installation is supported by 
GNU make utility. 

Installation of HEPGEN requires compiler appropriate for FORTRAN 77 standard, e.g. g77. In 
addition, installation requires CERNLIB library to be installed. After standard compilation the binary 
executable is created, hepgen . exe, which loads libraries dynamically. 



5.2 Input and output files 

HEPGEN use three input files - hepgen . data, fort . 9 and fort . 31. The text file hepgen . data 
is the configuration file of HEPGEN. It contains steering data, like for instance selection of simulated 
process, or used range of kinematic variables. All available options are described in Sec. 15.31 Exam- 



ple of hepgen . dat file is given in Appendix A The text file fort . 9 contains so-called random 



seed, i.e. integer number used for an initialization of a random number generator. The file is updated 
by the generator after finishing the job. The binary file f ort . 31, referred to as the beam file in the 
following, contains a set of beam particles. Each entry in the beam file corresponds to a single beam 
particle. An entry contains information about direction at a given space point, momentum and a type 
of beam particle which can be either true beam particle or beam halo particle. Simulation can be done 
either for true beam particles or for halo particles, or for both types (cf. Sec. 15.31) . The beam file is 
read by FORTRAN unformatted I/O stream mechanism. Instruction how to create a new beam file 
can be found in create Joe am . F file, which is enclosed to the HEPGEN source code. 

HEPGEN output files are fort . 11 and fort .4 0. The binary file fort . 11 contains output 
data. Access to this data is described in Sec. 15.41 The binary file fort . 4 contains histograms in 
HBOOK format, e.g. distributions of generated kinematics variables. 

The log of HEPGEN goes to the standard output. It contains for instance values of steering 
parameters and JETSET output of the first 30 events, which is useful for checks. 



5.3 Steering parameters 

All steering parameters are summarized in Tab. [TJ 
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Table 1: Description of steering parameters of HEPGEN. For each parameter the type (int/real) 
of argument (or arguments) is given and a list of possible values for those arguments which can 
take only specific, predefined values. 



Name 



Type of 
argument(s) 



Possible values Description 



NGEV int 
PROC int 

VMES int int 



DIFF 



int 



BMRD int 



BEAM real real int 
int int 



BPAR real real 
MACC real 

TPAR real real real 
real 



ALFA real 



a 




1 

2 
3 
4 
5 
7 





22 

211 

321 

-13 

-11 

211 




1 



1 

2 
3 



number of events to be generated 

type of process to be generated 

hard single production mode 

codes of generated particle and its decay mode 

7 (DVCS and BH) 

vr° — > 7 + 7 

p° -> 7T+ + 7T~ 

4> -> K + + K~ 
e+ + e~ 



{under development) 
{under development) 
{under development) 



J/* 
J/* 

p+ -> 7T+ + 7T° 

if set, events with nucleon diffractive dissociation are also 
generated 

diffractive dissociation off (only exclusive events) 
diffractive dissociation on (both exclusive and diffractive dis- 
sociation events) 
type of beam particles 

beam file not used, direction of the beam is parallel to z axis 
and momentum of the beam is given by the first argument of 
BEAM 

true beam only 
beam halo only 
true beam and beam halo 

definition of colliding particles, first and second arguments 
refer to the beam and the target energies, the third one is a 
PID of the beam particle, two last refer to the mass number 
and atomic number of the target, respectively 
lepton beam charge and polarisation 

azimuthal acceptance of scattered muons {this option is ob- 
solete) 

target parameters, the first argument is an average value of 
mass number of the target nuclei, second is a fraction of co- 
herent events to be generated in the sample, two last argu- 
ments are values of slopes of t' distributions for coherent and 
incoherent events, respectively 

dependence of cross section on the mass number of target 
nuclei, a is defined as a a = & P ■ A 0-1 , where a a is the cross 
section per nucleon for the reaction on the nucleus of mass 
number A, a p is the cross section for the reaction on the free 
proton 
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CUTL real real real minimum and maximum values of Bjorken scaling variable 

real real real (xsj), fraction of the lepton energy lost in LAB (y), nega- 

real real real tive four-momentum squared of the virtual photon (Q 2 ), total 

real energy squared in the virtual photon - nucleon system (W 2 ), 

energy of the virtual photon in LAB (u), energy of scattered 
beam particle and azimuthal angle of scattered beam particle 
respectively (most of those parameters became obsolete, the 
active ones are those related to Q 2 and v variables) 
TLIM real real minimum and maximum values of if variable, defined as if = 

\t\ — to, where t is four-momentum transfer to the target and 
to is minimal value of [t| allowed by the kinematics of the 
event 

REGG real real real bo b\ xq parametrization of slope b of ^distribution as a function of 

XBj, parametrization is given by b = bo + 2 • b\ log(xo /%Bj) 
(these parameters are used only for DVCS production) 



5.4 Access to output data 

The output data from HEPGEN generator are stored in fort . 1 1 file. The file can be used either by 
COMGEANT or directly by the user. In the first case, informations from the generator are available 
in PHAST, which is the framework for the COMPASS data analysis. Details how to access these 
informations in PHAST are given in Ref. |[T2ll . In the case of direct access the user should use the 
interface subroutine read_interf ace . F, which is enclosed to the HEPGEN source code. Infor- 
mation given below is relevant only when direct access to HEPGEN output data is needed (not needed 
for PHAST users). 

The output file f ort . 1 1 contains header record and event records. Each event record contains 
data, which are specific for corresponding event. The data, both from header and event records, are 
stored by read_interf ace . F in the following FORTRAN common blocks 

LEPTOU, USERVAR, LUJETS, 

where LEPTOU, USE R V AR , L U JE T S are common blocks , defined as 

REAL CUT (14) , PARL (30) , X, Y, W2 , Q2 , U 
INTEGER LST (40) 

COMMON /LEPTOU/ CUT, LST, PARL, X, Y, W2 , Q2 , U 

REAL USERVAR (2 0) 

COMMON /USERVAR/ USERVAR 

INTEGER* 4 N,K(4000, 5) 
REAL* 4 P(4000,5),V(4000,5) 
COMMON/LUJETS/N, K, P, V 

Description of the most useful variables is given in Tab. [21 Note, that the format of common block 
LUJETS and the values of elements of K array are consistent with those used in JETSET. Detailed 
description of JETSET can be found in Ref. lfT3l . 
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Table 2: Description of the most useful variables of HEPGEN output. 



Variable name 






Type 


Description 


USERVAR(3) 






real 


weight of event 

for px pi ii^ivp ntintnti nrnHiiPtinTi wpicrtit pnrrp'innTifl'i 

to|DVCS| 2 + |BH| 2 + Int 


TISFRVARCIfit 






1 v. el l 


for pypIiwivp nhotnn nrnrliiptinri wpioht for 1 TYVf", SI 1 ^ 

1U1 V_/ A 1 LI i> 1 V G U11ULU11 Ul UUU\^ L1U11 W^lgllL 1 V ^ 1 \M-J V \J kJ \ 


UOCl\V/\J\^l / ) 






real 


f~fw py r*] 1 1 ci \/p t~\ h 1 ti t~\ i*n (\ n i (\ n n/pi frl'it fir*T iT-^TTl^ 
1UI CACIUMVC Ullulull UIUU.UCL1U11 WClglll 1UI LDX1 


T KT( r )A\ 

l^O 1 ) 






inx 


A/PC ill IT Hi TTT*Qr*l"1 \ ' /"» Hi C C 1 Q ti (~\T\ \ ' t1 1 t1 / A | (li -i-f ov 

yCs v^/ Willi 11 VC lilaaUL-lallUll CVClll, 11U yVJ) 11 CA - 
p1l1<iivP PVPTlt 


X 






real 


value of Bjorken scaling variable (xsj) 


Y 
1 






real 


VdlUC Ul ItUlUll CllCLgy 1USI 111 L^r\D yy ) 


w? 

vv z. 






T*P5I 1 


vnliif* of total pnprov cmiarpH in trip virtual nhoton - 

VCllLlt Ul lULCll tllGlSiy >U Llill 111 Lilt VllLLltll U11ULU11 
niirOpon CVCtPTYl ^T/T/^^ 

nui^icuii sy&LCiii yvv ) 










volnp r\x npcrativp fni ir-TnnmpTi ti 1m cniidrpn r»f tnp vir- 

VCllLlC Ul llCgd-LlVC 1UU1 111U111C111L1111 dUUdlCU Ul L11C Vll 

tual nhoton (Cfi\ 

LLlal U11UIU11 ) 


TT 
U 






real 


valnp nf pnprcrv nf trip virtual nnntnti in T AR ( 7 i\ 

VtllLlC Ul CllClliy Ul L11C VllLUCll U11ULU11 111 LA 13 yis J 


N 






irrt 


nntrinpi* nf initial ititprtTiprli dtp 'i n rl rinsil ncirtir*lpc 

11U111UC1 Ul llllLldl, lllLCllllCUltlLC tlllLI lllldl Ual LltlCS 


K(T ^^ 


T - 1 


N 


VPS} 1 


QtatiiQ pnHp of 7^-th nartiplf* 

l>LuLL1l> LUUt Ul 1 111 UCULlL-lt 


K(7, 2) 


1= 1 


... ,N 


real 


PID of I-th particle, possible codes are summarized 
in Ref . HI 


K(7, 3) 


1= 1 


... ,N 


real 


index of parent particle of particle / 


K(J, 4) 


7= 1 


... ,N 


real 


index of first daughter of particle / 


K(/, 5) 


/= 1 


... ,N 


real 


index of last daughter of particle / 


P<7, 1) 


/= 1 


... ,N 


real 


momentum in x direction of particle / 


P(I, 2) 


/= 1 


... ,N 


real 


momentum in y direction of particle / 


P(J, 3) 


/= 1 


... ,N 


real 


momentum in z direction of particle / 


P(J, 4) 


/= 1 


... ,N 


real 


energy of particle / 


P(7, 5) 


/= 1 


... ,N 


real 


mass of particle / 
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Appendix A Example of hepgen. data file 



LIST 



* number of events 
NGEV 10 00 

* physics process 

PROC 5 

* produced particle and its decay mode 
VMES 

* diffractive dissociation switch 
DIFF 1 

* beam file 
BMRD 1 

* colliding particles 

BEAM 160.0 0.93827 -13 1 1 

* beam charge and polarisation 
BPAR 1.0 -1.0 

* acceptance of scattered muon 
MACC 0.050 

* target parameters 

TPAR 1.0 -0.1 52.2 5.0 

* dependence of cross sections on mass number A 
ALFA 0.7 5 

* limits for x B j, U, Q 2 , W 2 , v, E', (j) variables 

CUTL 0.0001 1.0 0.0 1.0 0.5 80.0 0.00 1000.0 

5.0 155.0 0.0 200.0 0.00 6.28318 

* limits for t' variable 
TLIM 0.005 1.0 

* parameterisaton ofxsj, t correlation 
REGG 4.94116 0.042 0.8 



END 
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